Quantum collapses and revivals of a matter wave in the dynamics of symmetry 
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Using the remarkable mathematical construct of Eugene Wigner to visualize quantum trajectories 
in phase space, quantum processes can be described in terms of a quasi-probability distribution 
analogous to the phase space probability distribution of the classical realm. In contrast to the 
incomplete glimpse of the wave function that is achievable in a single shot experiment, the Wigner 
distribution, accessible by quantum state tomography, reflects the full quantum state. We show that 
during the fundamental symmetry-breaking process of a generic quantum system - with a symmetry 
breaking field driving the quantum system far from equilibrium - the Wigner distribution evolves 
continuously with the system undergoing a sequence of revivals into the symmetry unbroken state, 
followed by collapses onto a quasi-classical state akin the one realised in infinite size systems. We 
show that generically this state is completely delocalised both in momentum and in real space. 

PACS numbers: 64.60.Ht, 05.30.-d,05.70.Ln 



I. INTRODUCTION 

Spontaneous symmetry breaking causes a macroscopic 
body under equilibrium condition to have less symme- 
try than its microscopic building blocks. Probably the 
phenomenon of superconductivity is the most spectacu- 
lar example of the symmetry breaking which a macro- 
scopic body spontaneously undergoes. Of course this is 
not the only one. Antiferromagnets, liquid crystals and 
other states of matter obey this rather general scheme of 
broken symmetries. The general idea is that when the 
number N of microscopic quantum constituents which, 
depending on the system, corresponds to the number of 
Cooper pairs, particles or spins, goes to infinity, the mat- 
ter undergoes a "phase" transition to a state in which 
the microscopic symmetries are violatecP^. In the con- 
text of quantum magnetisnP^, the macroscopic classical 
state has been described as a combination of "thin spec- 
trum" states emerging in the N —> oo limit because of 
the singular nature of the thermodynamic limit. The 
same description has been shown to apply also to the 
cases of quantum crystals, Bose-Einstein condensates and 
super conduct or s^HES The theory of spontaneous symme- 
try breaking explains the stability and rigidity of states 
which are not eigenstates of the underlying microscopic 
Hamiltonian but it makes no assertion on whether or 
how a symmetry broken groundstate can evolve out of 
the symmetric state in a closed quantum system. To in- 
vestigate this, one can perform a Gedankenexperiment 
in which a symmetry breaking perturbation is slowly 
switched on in an arbitrary large but finite system ini- 
tially prepared in a fully symmetric state. Using a partic- 
ular antiferromagnetic model system - the Lieb-Mattis 
model 6 - it has been recently shown that the corre- 
sponding quantum dynamics is dominated by highly non- 
adiabatic processes triggering the appearance of a sym- 
metric non-equilibrium state that recursively collapses at 



punctured times into a symmetry broken stated 

Here we shed light on this far-from-equilibrium 
symmetry-breaking process that is so ubiquitous in 
physics - in the formation of crystalline matter, atomic 
condensates, Josephson junction arrays and local pair- 
ing superconductors to name but a few - by introducing 
quantum state tomography. In complete analogy with 
medical diagnostics where three-dimensional images of 
the inner part of a body can be extracted from NMR or 
X-rays two-dimensional images obtained at different di- 
rections, quantum state tomography determines a quasi- 
probability distribution in phase space from only position 
(Q) or momentum (II) measurements^^. This quasi 
probability distribution has been introduced by Eugene 
Wigner in his phase space formulation of quantum me- 
chanics. For a pure quantum state, the Wigner distri- 
bution W(Q, II) is defined in terms of the position wave 
function *(Q) as W(Q, II) = tt" 1 **(Q - S)9(Q + 

S)e~ 2,lUS dS (in ft = 1 units) and retains the marginal 
probability distributions^ 

J W(Q,U)dU= |4>(Q)| 2 , 

J W(Q,U)dQ = \^>(U)\ 2 . 

By detecting the position of many objects prepared in 
the same quantum state yields the spatial distribution 
of the wavefunction |^(Q)| 2 as a space-like shadow of 
the Wigner distribution. This, in turns, allows for a to- 
mographic reconstruction of W(Q, II) once various shad- 
ows at different directions in phase space have been ob- 
served. Therefore, although humans are chained in a 
Plato's quantum cave and are entitled to see only shad- 
ows of a quantum object, they can access full quantum 
snapshots once these shadows are properly rotated^. 
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We unravel these snapshots in the dynamics of sym- 
metry breaking by using the paradigmatic example of an 
harmonic crystal to show that the quantum dynamics is 
generically characterized by the appearance of revivals 
of the symmetric ground state wavefunction followed by 
collapses towards a quasi-classical state akin the symme- 
try broken groundstate of infinite size system. In this 
quasi-classical state, however, the matter wave has max- 
imum uncertainty both in total position - precisely as in 
the symmetric translational invariant ground state - and 
in total momentum. The exceptions are punctured times 
which render a Dirac comb of symmetry broken states^. 
Interestingly we find this sequence of collapses/revivals of 
the ground state wavefunction to occur on a characteris- 
tic timescale set by Zurek's equation of non-equilibrium 
quantum phase transition^. 



II. 



SPONTANEOUS SYMMETRY BREAKING 
AND THE THIN SPECTRUM. 



The Nambu-Goldstone theorem^ 20 * 21 * guarantees the ex- 
istence of low-energy gapless collective excitations in sys- 
tems with spontaneously broken continuous symmetries. 
The low-energy Hamiltonian for these normal modes, 
which depending on the particular system at hand corre- 
spond to phonons, spin waves or Bogoliubov's excitations 
can be always recast in the form 



(i) 



where c^k indicate the frequencies of the Goldstone mode 
excitations vanishing in the long- wavelength limit. If 
there is a spontaneously broken symmetry, the motion 
along the continuous symmetry axis characterizing the 
quantum mechanics of the macroscopic body as whole 
cannot be of the form 12 b^b but it will be either given 
by n 2 /(27) (as in crystals and Josephson junction ar- 
rays) or Q 2 1 (21) (e.g. antiferromagnets and local pairing 
superconductors)^^ where Q and II are the coordinate 
and the conjugate momentum operator along the sym- 
metry axis whereas the parameter / depends on the to- 
tal number of microscopic quantum constituents and di- 
verges in the thermodynamic limit where TV — )• oo. This 
is explicitly manifest by considering the specific example 
of an harmonic crystal with Hamiltonian 



(2) 



where the index j labels the N atoms in the lat- 
tice which have mass m, momentum pj, and posi- 
tion Xj. The harmonic potential among neighboring 
atoms is parametrized by the constant n. Let us now 
define the bosonic annihilation and creation operator 
bj = [Cxj + ipj/C] /V2h and b\ = [C xj - ipj/C] /V2h 



where C = (2m^) 1//4 . In momentum space, the Hamil- 
tonian Eq.Q reduces to 



n = [Mh + ^r(bW_ k + b k b. k ) + 1 



where = 2 — cos(/ca), B k = — cos(ka) and a is the lat- 
tice constant. This Hamiltonian is not diagonal since the 
terms b\b^_ k and bkb-k create and annihilate two bosons 
at the same time. One can get rid of them by perform- 
ing a canonical Bogoliubov transformation. However, the 
parameters in the Bogoliubov transformation diverge as 
k —> and thus the canonical transformation is not well 
defined^. This implies that one should investigate sepa- 
rately the k = component 



H«-<») ! 



This part of the Hamiltonian describes the fact that the 
quantum crystal carries a kinetic energy associated with 
the combined mass of all N atoms. Going back to real 
space it reads 



n 2 



k=0 



2m7V' 



(3) 



where n = ^2jPj is the total momentum of the entire 
system and we left out a negligible constant. The ground 
state of this Hamiltonian has total zero momentum: it 
has no uncertainty in total momentum and maximum 
uncertainty in total position thereby implying that the 
translational symmetry is unbroken. At finite TV, the ex- 
citations over the ground state respecting the symmetry 
give rise to a tower of ultra- l ow en ergy states that form 
the so-called "thin spectrum"! 10 ! 11 ! . It is "thin" because 
it contains states that are of such low energy that their 
contribution to thermodynamic quantities vanish in the 
thermodynamic limit. Nevertheless, when N — »• oo, the 
thin spectrum excitations collapse to form a degenerate 
continuum of states. Within this continuum even a van- 
ishingly small symmetry breaking perturbation is enough 
to couple different thin spectrum states thereby stabiliz- 
ing the symmetry broken ground state. To show this, 
let us take into account the effect of a symmetry break- 
ing perturbation - a pinning potential for the individ- 
ual atoms- rendering a symmetry breaking Hamiltonian 
Hsb — —B cos (^TTXj) j (27r) 2 . For small deviations of 
the atoms from their mean positions, the zero momentum 
term in lowest order is giverP^ by %sb = Bx 2 ot / '(27V) 
where x tot = ^ Xj. This, in turn, implies that the col- 
lective behaviour of the harmonic crystal as a whole is 
governed by the harmonic-oscillator Hamiltonian 



2N 



BNQ 2 



(4) 



where for simplicity we have set m = 1 and we intro- 
duced the centre of mass coordinate Q = x to t/N sat- 
isfying the canonical commutation relation [Q,II] = ih. 
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FIG. 1: (Color online), (a) Sketch of the behavior of the 
relaxation time (continuous line) as compared to the typical 
time on which the Hamiltonian is changed t determining the 
freeze-out time of Zurek's equation, (b) The corresponding 
regimes for the dynamics of symmetry breaking in the t — to 
plane. Times have been measured in units of the freeze-out 
time t. The inset shows the setup of the symmetry breaking 
field. 



The quantum of energy of this Hamiltonian oc yB and 
the excitations over the ground state realize a "dual" 
thin spectrum. The ground state wavefunction corre- 
sponds to a Gaussian wave packet for the collective co- 
ordinate Q of the form ^o(Q) e_< ^ ^ 2L with width 
/ r~\~ 1/2 

L oc IN v B ) . For a vanishing symmetry break- 
ing field and a finite number of atoms, we find that the 
ground state wavefunction obviously collapses onto the 
symmetric ground state of the microscopic Hamiltonian. 
However, by taking first the thermodynamic limit the 
centre of mass position becomes completely localised even 
if at the end the symmetry breaking field is sent to zero. 
Therefore one finds that the system is in a stable state 
which is not an eigenstate of the underlying microscopic 
Hamiltonian - the system is inferred to spontaneously 
break the symmetry. Strictly speaking, only truly in- 
finite size systems are allowed to spontaneously break 
the symmetry. A large, but not infinitely large crystal 
requires a finite symmetry breaking field to stabilize a 
symmetry broken state over the symmetric ground state 
of the microscopic Hamiltonian. 



III. QUANTUM DYNAMICS OF SYMMETRY 
BREAKING 

Let us then consider such a large but finite system with 
a pinning potential whose strength is switched on linearly 
in time as B(t) = 5t with ramp rate S. At initial time 
to we consider a finite symmetry breaking perturbation 
B (to) = Bq [ c -f- inset of Fig. [ljb)] and the wave func- 
tion of the system in the corresponding ground state. We 
stress that the choice of a finite symmetry breaking per- 
turbation at initial time is essential in order to have a 
cutoff guaranteeing the continuity of the wave function 
basis. Lateron we will consider the limit Bo — >• corre- 
sponding for finite N to an initial completely symmet- 
ric ground state with gapless dual thin spectrum excita- 
tions. The inclusion of a cutoff Bo renders two distinct 
regimes of the quantum dynamics. Whenever the char- 

1/2 

acteristic relaxation time r = (St) is much smaller 
than the typical timescale t on which the Hamiltonian 
is changed, the system is able to react to the changing 
Hamiltonian thereby rendering an adiabatic passage. In 
the opposite regime r ^> t, reflexes of the system are 
so deteriorated that the state can be considered effec- 
tively frozen and the dynamics is impulse-like. Clearly, 
the crux of this story is the instant where the dynam- 
ics changes from impulse to adiabatic -it is determined 
by Zurek's equation^ r(t) = t and defines the so-called 
freeze-out time i oc c) -1 / 3 [c.f. Fig. [lja)]. By considering 
an initial time to ^> t, the entire evolution will be thus 
nearly adiabatic. In this case, fluctuations of the cen- 
tre of mass coordinate decreases continuously in time as 

-l 



AQ 2 = (Q 2 ) - (Q) 2 oc 



However, for the dy- 



Ny/St 

namics of symmetry breaking to be adiabatic the ramp 

3/2 

rate 5 is seen to be bounded by S < B f . Henceforth 
for a vanishing ramp rate S at finite values of £?o, we re- 
cover a quasi- adiabatic time evolution. But taking the 
Bo —> limit at finite ramp rate 5, the dynamics will 
start in the impulse regime even if at the end the ramp 
rate is sent to zero - the adiabatic limit can never be 
reached for a sufficiently small Bq. This is in agreement 
with the recent finding 2 ^ that adiabatic processes in low- 
dimensional systems with broken continuous symmetries 
are absent. 



In the same spirit of Ref. [T4] to analyse the dynamics of 
symmetry breaking in the strongly non- adiabatic regime, 
we first make use of the adiabatic-impulse (AI) approx- 
imation which underlies the Kibble-Zurek theor y 23 " 25 ! of 
non-equilibrium phase transition. In the AI scheme^, 
the initial state is considered effectively frozen in the 
impulse-like regime to < t < t changing only by a 
trivial overall phase factor. At freeze-out time, the 
system therefore reaches a state that is a superposi- 
tion of dual thin spectrum excitations |^o(Q^o)) = 
c n \ty n (Q,t)) where the coefficients c n are non zero 
only for even values of the quantum number n. The 




FIG. 2: (Color online) (a) Time evolution of (AQ 2 ) -1 for different values of the initial time to- The continuous thick line is 
for to = 10 _1 , the dashed line for to = 10 -2 and the dotted line corresponds to to = 10 -3 . By decreasing the initial time to, 
a Dirac comb of symmetry-broken states is approached, (b) Same for the inverse of the fluctuations of the total momentum. 
Also in this case a Dirac comb is approached, (c) Same of (a) for the time evolution of the total uncertainty of the matter 
wave. Apart from punctured times, it is seen to diverge in the to — »• limit. 



evolution at t > t can be considered to be adiabatic 
and therefore the dynamics of the wavefunction is gov- 
erned by |tf(Q,*)) = En c nl*n(Q,«))e- <n »W where 
we have defined the dynamical phase factor Q n (t) = 

2/3 x (t/if /2 - 1 (n + 1/2) . Within the AI approxi- 
mation, we can obtain the time-evolution of the symme- 
try breaking order parameter - defined by the inverse of 
the fluctuations of the centre of mass coordinate- tak- 
ing explicitly into account quantum phase interference 
effects. This is unlike for instance the case of the Lieb- 
Mattis model where in the AI scheme the time evolution 
of the staggered magnetisation can be computed only by 
neglecting interference effects^. We find 



of the total momentum of the entire crystal and find 



[AQ 2 (t)]" 1 = 2^ 




sin 



(5) 

the behaviour of which is shown in Fig. ^[a). When de- 
creasing the initial time to we observe that the behaviour 
of the order parameter calculated above corresponds pre- 
cisely to a Dirac comb of symmetry broken states in per- 
fect agreement with the case of antiferromagnets 14 . Be- 
sides this, we find that at the punctured times where 
the fluctuation of the centre of mass coordinate van- 
ishes the dynamical phases of the excited dual thin 
spectrum excitations have tt shifts, i.e. for t = t& = 
[3ft7r/2 + 37r/4 + l] 2 ^ 3 with k integer. This shows that in- 
dependent of the actual strength of the symmetry break- 
ing perturbation, destructive quantum phase interference 
leads to an instantaneous breaking of the translational 
symmetry. It is also manifested by the fact that a direct 
computation shows at these instants Q\Sfr(Q,tk)) = 
and therefore the harmonic crystal is completely localised 
in the centre of the potential well. 

To further show the nature of the non-equilibrium state 
realised in the remaining time evolution, we have deter- 
mined the time-evolution of the inverse of the fluctuations 



[An 2 (t)]- 




i 



i 



sin' 



(6) 

the behaviour of which is shown in Fig. |2|b). A Dirac 
comb is also the result. The instants where the system 
is an eigenstate of the total momentum II - which cor- 
respond to revivals of the initial completely delocalised 
symmetric state even in the presence of a sizable pinning 
potential - are realised for t = t/~ = [3^7r/2 + l] 2 / 3 with n 
integer in which case quantum phase interference effects 
are absent. This is again in line with the dynamics of the 
Lieb-Mattis model and henceforth guarantees the univer- 
sality of the dynamical symmetry breaking phenomenon 
independent of the specific microscopic model taken into 
account. 

Apart form the punctured times where the matter wave 
has either no uncertainty in total position or no uncer- 
tainty in total momentum, a strongly non-equilibrium 
state is realised. This state is rather interesting as it re- 
tains a perfect delocalisation of both the centre of mass 
position and its corresponding momentum. It has an in- 
finite uncertainty, i.e. An 2 AQ 2 — » oo as it is shown in 
Fig. [2b. 



IV. QUANTUM DYNAMICS IN PHASE SPACE 

To unravel the origin of this non-equilibrium state we 
analyse the dynamics of symmetry breaking with quan- 
tum state tomography. The time-dependent Hamiltonian 



U(t) 



2N 



N 



StQ 2 , 



represents a simple example of generalized time- 
dependent harmonic oscillator whose exact quantum the- 
ory has been extensively studied in the literatur e 27 " 39 !. In 
particular, within the Feynman path integral approach it 
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FIG. 3: (Color online) Density plots of the Wigner function in 
phase space (Q, II). The centre of mass coordinate Q and the 
total momentum II have been rescaled by iV _1//2 and N 1 ^ 2 re- 
spectively in order to absorb the N dependence of the Wigner 
function. By decreasing the initial time to the time evolution 
changes from an adiabatic shear motion, to a strongly non- 
equilibrium rotative motion. 



has been showrP^ that the spectral decomposition of the 
propagator Q(Qb,t b \Q a ,t a ) = J2 n ^n(Qa,t a )^n(Qb,t b ) 
is defined by a complete set of wavefunctions of the form 



*n(Q,*) = 



2 n n\ 



Re [u(t)] 



1/4 



^Re [u>(t)] Q 



x e 



i(n+|)0(t) 



(7) 



where H n are the Hermite polynomials and Re [w(t)] > 
in order to guarantee square integrability. The quantal 
phase <p(t) and the complex parameter uo(t) are uniquely 
determined by solving the classical Euler-Lagrange equa- 
tion of motion. Different sets of wave function of the 
form Eq. [7] correspond to take different pairs of linearly 
independent solutions to the classical equation of motion. 
This enables to choose two particular solutions guaran- 
teeing that at the initial time ^o(Q?^o) corresponds to 
the initial static Gaussian wavepacket and implies that 
the wavefunction at all times remains an n = state of 
the form of Eq. [F| 

The exact solution of the time-dependent Schrodinger 
equation allows us to determine the time evolution of the 
Wigner function [see Fig{3] given by 



W(Q,U,t) = ± e - Re[Q{t)]Q2 x e -(n+/m[fi(t)]Q)Vfle["(*)]. 

7T 

(8) 

By considering a nearly-adiabatic process - an initial 
time to » t - one finds that the Wigner function shears 
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FIG. 4: (Color online) The time evolution of the angle 6 char- 
acterising the rotation in phase space of the Wigner function. 
By increasing the number of microscopic constituents, a step- 
like behaviour is realised. 



in time in agreement with the time evolution of a Gaus- 
sian wavepacket reacting adiabatically to the time change 
of the harmonic oscillator angular frequency. On the con- 
trary, in the out-of-equilibrium regime, i.e. for to < t, the 
Wigner distribution shears and rotates in phase space as 
it follows from the fact that the initial width of the Gaus- 
sian wavepacket acquires a non- negligible imaginary part. 
Finally, in the to — >• limit instead, any shear is absent 
and the motion simply corresponds to a rigid rotation in 
phase space. For a completely symmetric initial state, the 
initial Wigner function corresponds to an infinite line in 
phase space S(Q) and a rigid rotation is the only motion 
preserving this one-dimensional character. As a result, 
we find in the to —> limit the time evolution of the 
Wigner function as 



W(Q, n, t) ~ £(cos 0(t) Q + sin 0(t) Ft), 



(9) 



where the angle 0(t) ~ tan -1 Im[Q(t)]. By increasing 
the number of microscopic quantum constituents N, the 
time dependence of the angle 0(t) approaches a step-like 
behavior [c.f. Fig. [4] thereby implying that the quan- 
tum dynamics of symmetry breaking in a macroscopic 
body is characterized by revivals of the initial symmetric 
state and collapses of the initial quantum state towards 
a "quasi-classical" state - tomographically indistinguish- 
able from the symmetry broken state of infinite size sys- 
tems but, as we have shown above, completely delocalised 
both in momentum and in real space. 



V. CONCLUSIONS 

In conclusion, by considering the paradigmatic exam- 
ple of an harmonic oscillator, we have shown that in the 
dynamical realm symmetry breaking is characterized by 
far-from-equilibrium processes. No matter how slowly 
a symmetry breaking perturbation is driven, the adia- 
batic limit can never be reached in a macroscopic body. 
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By means of quantum state tomography, we have shown 
that nevertheless the evolution of symmetry breaking cor- 
responds to a continuous, rigid rotation of the Wigner 
distribution. This rotation yields at the same time a 
sequence of step-like revivals of the symmetric state fol- 
lowed by collapses onto a symmetry broken groundstate 
akin to the one realized in infinite size systems but with 
maximum uncertainty both in total position and in total 



momentum. 
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